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We investigate if known extrinsic and intrinsic factors fully account for the complex features observed in 
recordings of human activity as measured from forearm motion in subjects undergoing their regular daily rou- 
tine. We demonstrate that the apparently random forearm motion possesses previously unrecognized dynamic 
patterns characterized by fractal and nonlinear dynamics. These patterns are unaffected by changes in the aver- 
age activity level, and persist when the same subjects undergo time-isolation laboratory experiments designed 
to account for the circadian phase and to control the known extrinsic factors. We attribute these patterns to a 
novel intrinsic multi-scale dynamic regulation of human activity. 

PACS numbers: 87.19.st, 87.80.-y, 87.90.+y, 89.20.-a 



Control of human activity is complex, being influenced by 
many factors both extrinsic (work, recreation, reactions to 
unforeseen random events) and intrinsic (the circadian pace- 
maker that influences our sleep/wake cycle [1] and ultradian 
oscillators with shorter time scales [2]). The extrinsic fac- 
tors may account for the apparently random fluctuations in hu- 
man motion observed over short time scales while the intrinsic 
rhythms may account for the underlying regularity in average 
activity level over longer periods of up to 24 h. Further, hu- 
man activity correlates with important physiological functions 
including whole body oxygen consumption and heart rate. 

Actiwatch devices [3] are traditionally used to demarcate 
sleep versus wakefulness based on average activity levels, or 
to observe the mean pattern of activity as it changes across the 
day and night according to disease state (Fig. 1). Tradition- 
ally activity fluctuations are considered as random noise and 
have been ignored. We hypothesize that there are systematic 
patterns in the activity fluctuations that may be independent 
of known extrinsic and intrinsic factors. 

To test our hypotheses, we evaluate the structure of human 
activity during wakefulness, using; (i) probability distribution 
analysis; (ii) power spectrum analysis, and (iii) fractal scaling 
and nonlinear analysis. To elucidate the presence of an in- 
trinsic activity control center independent of known circadian, 
ultradian, scheduled and random factors, we apply 3 comple- 
mentary protocols. 

• (A) Daily routine protocol: We record activity data 
throughout two consecutive weeks in 16 healthy ambulatory 
domiciliary subjects (8 males, 8 females, 19-44 years, mean 
27 years) performing their routine daily activities. The only 
imposed constraints are that subjects go to bed and arise at the 
same time each day (8 h sleep opportunity) and that they are 
not permitted to have daytime naps (Fig. 1). 

• (B) Constant routine protocol: To assess intrinsic activity 
controllers (i.e. circadian or other neural centers) independent 
of scheduled and random external influences activity record- 
ings are made in the laboratory throughout 38 h of constant 
posture (semi-recumbent), wakefulness, environment (21°C, 



dim light [< 8 lux]), dietary intake and scheduled events [4]. 
This protocol is performed in a subset of subjects (7 males, 4 
females) that participated in the daily routine protocol. These 
highly controlled and constant experimental conditions result 
in reduced average and variance of activity levels. 

• (C) Forced desynchrony protocol: To test for the pres- 
ence of heretofore unidentified intrinsic activity control cen- 
ters, independent of known activity regulators (circadian pace- 
maker), whilst accounting for scheduled and random external 
influences, we employ the vaUdated Forced desynchrony (FD) 
protocol [1]. Six (4 male, 2 female) of the 16 subjects that par- 
ticipated in the daily routine protocol completed the FD limb 
of the study. For eight days subjects remain in constant dim 
light (to avoid "resetting" the body clock). Sleep periods are 
delayed by 4 h every day, such that subjects live on recurring 
28 h "days", while all scheduled activities become desynchro- 
nized from the endogenous circadian pacemaker. Thus, as 
measurements occur across all phases of the circadian clock, 
the effect of intrinsic circadian influences can be removed [1]. 
Average activity level and activity variance are also signifi- 
cantly reduced due to laboratory-imposed restrictions on the 
subjects activity (Fig. 1). 

When the same subject is studied in different protocols, we 
find large differences in the probability distributions (Fig. 2). 
For example, during wakefulness greater values of activity 
occur most frequently during the daily routine, intermediate 
activity values occur during the forced desynchrony, and the 
highest frequency of low activity values is seen during the 
constant routine (Fig. 2a). Indeed, the largest activity val- 
ues encountered during the constant routine protocol are ap- 
proximately two orders of magnitude less frequent than sim- 
ilar activity values encountered in the daily routine protocol. 
We also find major differences between individuals in the dis- 
tribution of activity values during the daily routine protocol 
(Fig. 2b). Such differences are expected given their differ- 
ent daily schedules, environments, and reactions to random 
events. However, by appropriately rescaling the distributions 
of activity values on both axes to account for differences in 
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FIG. 1 : Independent contributors to the complex dynmaics of hu- 
man activity, depicted at the top of the figure, include: Q) reaction to 
extrinsic random events, @ scheduled activities and, (3) intrinsic fac- 
tors, notably the endogenous circadian pacemaker which influences 
the sleep/wake cycle. Our findings of scale-invariant activity patterns 
(Figs. 2, 3) indicate a heretofore-unidentified intrinsic multi-scale 
control of human activity (J), which is independent of other extrinsic 
and intrinsic factors such as ®, @, and (3). The second panel illus- 
trates an actual one-week recording of human activity [3] during the 
daily routine protocol. Data structure highlights a 24-h sleep/wake 
periodic change in the mean activity — lowest during sleep (filled 
bars). The third panel, expanding a I6-h section of wakefulness, also 
shows patches of high and low average activity levels with apparent 
erratic fluctuations at various time scales. The bottom left panel is 
an activity recording from the same subject during the constant rou- 
tine protocol with much lower average activity values compared to 
daily routine. The clear 2-h cycle is a result of scheduled labora- 
tory events. The bottom right panel shows activity levels in the same 
subject during the forced desynchrony protocol. Of note here is the 
28-h sleep/wake cycle as opposed to the 24-h rhythm in activity data 
during the daily routine. 

average activity level and standard deviation [5], we find a 
remarkable similarity in the shapes of the probability distri- 
butions for all three protocols (Figs. 2c, 2e), and for all indi- 
viduals when in the same protocol (Fig. 2d). The existence of 
a universal form of the probability distribution, independent 
of activity level in all individuals and in all protocols, sug- 
gests that a common underlying mechanism may account for 
the overall distribution of activity. 

This probability distribution when plotted on a log-log scale 
reveals different characteristics above and below a distinct 
crossover point (Fig. 2e). At scales above the crossover ac- 
tivity level there is pronounced non-Gaussian tail (Fig. 2e). 
This tail on the log-log plot represents a power-law form, in- 
dicating an intrinsic self-similar structure for a range of ac- 



FIG. 2: Probability distributions of activity values during wakeful- 
ness, (a) Probability distributions for an individual subject during 14 
consecutive days of daily routine, 38 h of constant routine and 8 days 
of the forced desynchrony protocol, (b) Probability distributions for 
16 subjects during the daily routine protocol indicate large difference 
between individuals, (c) Same probability distributions as in (a), after 
appropriately rescaling both axes. Data points for all three protocols 
collapse onto a single curve, (d) Same probability distributions as in 
(b), after rescaling. (e) Group average of the rescaled distributions 
during all three protocols. Again all distributions collapse onto a sin- 
gle curve, (f) Group average of all individual distributions rescaled 
as in (e) obtained for varied time windows during the daily routine. 



tivity values. Moreover, we find that the observed shape of 
the rescaled probability distribution remains unchanged when 
the data series are reanalyzed using a variety of observation 
windows ranging from 15 s to 6 min (Fig. 2f). This stability 
of the probability distribution over a range of time scales in- 
dicates that the underlying dynamic mechanisms controlling 
the activity have similar statistical properties on different time 
scales. Statistical self-similarity is a defining characteristic of 
fractal objects and is reminiscent of a wide class of physical 
systems with universal scaling properties. Our finding of a 
universal form of the probability distribution raises the possi- 
bility of an intrinsic mechanism that influences activity values 
in a self-similar "fractal" manner, that is unrelated to the indi- 
vidual's daily and weekly schedules, reactions to the environ- 
ment, the average level of activity, the phase of the circadian 
pacemaker, and the time scale of observation. 

We next perform power spectral analyses for all three proto- 
cols to determine whether there exist any systematic intrinsic 
ultradian rhythms of activity with periods of less than 24 h du- 
ration [2, 6]. The data for each individual exhibit occasional 
peaks in the daily routine protocol for periods ranging from 30 
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FIG. 3: Power spectrum and correlation analyses of activity data 
during wakefulness, (a) Group average power spectral densities for 
all three protocols. For better clarity and to avoid overlap, curves 
are vertically offset. To present graphs on a common x-axis, power 
spectra are shown with decreasing frequency from left to right. The 
spectral density peaks for the simulated scheduled activity data rep- 
resenting controlled scheduled events during the protocol (bottom 
curve) match the peaks observed in the original human activity data 
recorded during the forced desynchrony protocol. This analysis sug- 
gests that the observed peaks in the power spectrum cannot be at- 
tributed to endogenous ultradian rhythms, (b) DFA scaling of ac- 
tivity for an individual during wakefulness throughout two separate 
weeks of daily routine, 38 h of constant routine, and 8 days of forced 
desynchrony protocols, (c) DFA scaling of the magnitude series of 
activity increments for the same signals as in (a). A scaling expo- 
nent Qfmag ~ 0.8 of similar value is observed for all three proto- 
cols, (d) Scaling exponents a and Qmag (left scale), and average 
activity levels (right scale) for all 16 subjects obtained from a 14-day 
daily routine protocol. Although the average activity level between 
subjects changes considerably (from 0.2 to 0.5), both scaling expo- 
nents are consistent for all subjects, exhibiting a group average of 
a = 0.92 ± 0.05 and a^ag = 0.77 ± 0.05. 



min to 4 h. However, we find no systematic ultradian rhythms 
within individuals from week to week, and no systematic ul- 
tradian rhythms in the group average for the daily routine pro- 
tocol (Fig. 3a). The only systematic rhythms that are ostensi- 
bly in the ultradian range which emerge in the group data are 
at 4 h during the forced desynchrony protocol (with harmonics 
at 2 h and 80 min) and at 2 h during the constant routine pro- 
tocol (with harmonics at 1 h and 30 min) (Fig. 1 and Fig. 3a). 
These peaks are caused by the controlled scheduled activities 
in the laboratory and are extrinsic to the body as they also oc- 
cur in simulated scheduled activity data that assumes specific 
activity values for each scheduled behavior imposed through- 
out the laboratory protocols (Fig. 3a). Thus, we find no evi- 
dence of systematic intrinsic ultradian rhythms in our data. 

To provide further insight into the dynamic control of ac- 
tivity, we next examine the temporal organization in the fluc- 
tuations in activity values. We perform detrended fluctuation 
analysis (DFA) which quantifies correlations in the activity 
fluctuations after accounting for nonstationarity in the data by 
subtracting underlying polynomial trends. The DFA method 
quantifies the root mean square fluctuations, F{n), of a sig- 
nal at different time scales n. Power-law functional form, 
F{n) ~ n", indicates self-similarity (fractal scaling). The 
parameter a, called the scaling exponent, quantifies the corre- 
lation properties in the signal: if a=0.5, there is no correlation 
(random noise); if a < 0.5, the signal is anticorrelated, where 
large activity values are more likely to be followed by small 
activity values; if a > 0.5, there are positive correlations, 
where large activity values are more likely to be followed by 
large activity values (and vice versa for small activity values). 

Figure 3b shows that F{n) for a typical subject during 
wakefulness exhibits a power-law form over time scales from 
Ri 1 min to w 4 h. We find that the scaling exponent a is 
virtually identical for records obtained during the first week 
of daily routine {a = 0.92 ± 0.04, mean ± standard deviation 
among subjects), the second week {a ~ 0.92 ± 0.06) of the 
daily routine, the constant routine protocol {a = 0.88 ± 0.05), 
and the forced desynchrony protocol {a = 0.92 ± 0.03). The 
value of a « 0.9 for all protocols and all individuals in- 
dicates that activity fluctuations are characterized by strong 
long-range positive correlations. Furthermore, we find that 
this scaling behavior is not caused by the scheduled activi- 
ties because simulated scheduled activity data that are gener- 
ated by assigning a specific activity value for each scheduled 
event throughout the laboratory protocols yields an exponent 
of a = 1.5 (Fig. 3b), which represents random-walk type be- 
havior These results suggest that the activity fluctuations are 
not a consequence of random events (in which case a would 
be 0.5) or scheduled events, but rather relate to an underlying 
mechanism of activity control with stable fractal-like features 
over a wide range of time scales from minutes to hours. Since 
mean activity levels and the amplitude of the fluctuations are 
greatly reduced in the laboratory during the constant routine 
and forced desynchrony protocols (Fig. 1), we obtain smaller 
values of F{n) (downward shift of the lines in Fig. 3b). How- 
ever there is no change in the scaling exponent a. Similarly, 
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the scaling exponents for the daily routine protocol are inde- 
pendent of the average activity levels of the different subjects 
(Fig. 3d), and of the mean activity level on different days of 
the week, indicating that this newly-found scaling pattern of 
activity appears to be an intrinsic feature. 

To test for the presence of nonlinear properties of the data, 
we analyze the "magnitude series" formed by taking the ab- 
solute values of the increments between consecutive activity 
values [7]. Again, from detrended fluctuation analysis of this 
series, we find practically identical scaling exponents, amag, 
for all three protocols, despite large differences in mean ac- 
tivity levels between protocols (Fig 3c). Moreover, all indi- 
viduals have very similar values of the scaling exponent ctmag 
(Fig. 3d), which are not systematically changed by the proto- 
col. For the group, during the first week of daily routine, we 
find ttniag = 0.78 ± 0.06 (mean ± standard deviation among 
subjects), during the second week amag = 0.76 ± 0.05, dur- 
ing the constant routine protocol amag = 0.82 ± 0.05, and 
during the forced desynchrony protocol a,nag = 0.80 ± 0.04. 
Since Q!,„ag ~ 0.8(> 0.5), there are positive long-range cor- 
relations in the magnitude series of activity increments, indi- 
cating the existence of nonlinear properties related to Fourier 
phase interactions (Fig. 3c) [7, 8]. To confirm that the ob- 
served positive correlations in the magnitude series indeed 
represent nonlinear features in the activity data, we do the fol- 
lowing test: we generate a surrogate time series by performing 
a Fourier transform on the activity recording from the same 
subject during daily routine as in Fig. 3b, preserving the am- 
plitudes of the Fourier transform but randomizing the phases, 
and then performing an inverse Fourier transform. This proce- 
dure eliminates nonlinearities, preserving only the linear fea- 
tures of the original activity recording such as the power spec- 
trum and correlations. Thus, the new surrogate signal has the 
same scaling behavior with a = 0.93 (Fig. 3b) as the original 
activity recording; however, it exhibits uncorrected behavior 
for the magnitude series (amag = 0.5) (Fig. 3c). Our results 
show that the human activity data contains important phase 
correlations which are canceled in the surrogate signal by the 
randomization of the Fourier phases, and that these correla- 
tions do not exist in the simulated scheduled activity. Further, 
our tests indicate that these nonlinear features are related to 
Fourier phase interactions, suggesting an intrinsic nonlinear 
mechanism [8]. The similar value of amag for all three proto- 
cols and all individuals, which is different from amag ~ 0.5 
obtained for the simulated scheduled activity and for the phase 
randomized data, confirms that the intrinsic dynamics pos- 
sess nonlinear features that are independent of the daily and 
weekly schedules, reaction to the environment, the average 
level of activity, and the phase of the circadian pacemaker 

The consistency of our results among individuals, and for 
different protocols, suggests that there exist previously unrec- 



ognized complex dynamic patterns of human activity that are 
unrelated to extrinsic factors or to the average level of activity. 
We also showed these patterns to be independent of known 
intrinsic factors related to the circadian and to any ultradian 
rhythms. Notably, (i) these patterns are unchanged when ob- 
tained at different phases of the circadian pacemaker; (ii) we 
do not observe systematic intrinsic ultradian rhythms in activ- 
ity among subjects in the daily routine experiment; (iii) im- 
posing strong extrinsic ultradian rhythms at 4 h and 2 h in the 
laboratory protocols did not change the fractal scaling expo- 
nents a or amag or the form of the probability distribution; 
and (iv) we find consistent results over a wide range of time 
scales. Together, these findings strongly suggest that our re- 
sults are unlikely to be a reflection of the basic rest activity 
cycles or ultradian rhythms. We attribute these novel scale- 
invariant patterns to a robust intrinsic multi-scale mechanism 
of regulation (Fig. 1). This regulatory mechanism presents a 
new challenge to understand nonlinear control of human mo- 
tor activity and pathways of interaction with other physiologic 
dynamics such as heart rate, gait [9], finger tapping [10], hu- 
man sway and muscle fluctuations [11]. 



[1] C. A. Czeisler etal. Science 284, 2177 (1999). 
[2] N. Kleitman, Sleep 5, 311 (1982). 

[3] The subject wears a wristwatch-sized Acf/wflfc/i recorder (Mini- 
Mitter Co., OR, USA) that unobtrusively measures changes in 
forearm acceleration in any plane (sensitive to O.Olg, where g is 
the acceleration due to gravity). See, e.g., G. Jean-Louis et at.. 
Physiology & Behavior 70, 49 (2000), and references therein. 
Each data point recorded in the device's internal memory repre- 
sents the value of changes in acceleration sampled at 32 Hz and 
integrated over a 15-second epoch length. Recordings are made 
for up to two weeks, yielding approximately 10^ data points for 
each subject. Inhomogeneity of recording sensititiy across the 
range of activities is accounted for in the analyses. 

[4] C. A. Czeisler et al. Science 244, 1328 (1989); E. N. Brown 
and C. A. Czeisler, J. Biol. Rhythms 7, 177 (1992). 

[5] To test if the individual probability density curves follow a com- 
mon functional form, we rescale the curves while preserving 
the normalization to unit area. We divide the activity values 
by a constant, Aq, and multiply the probability density func- 
tion by the same constant, where Aq is the activity value be- 
fore rescalingof each individual curve for which the cumulative 
probability (i.e., the area under the density function curve) is 
60%. 

[6] B. D. Lewis et al.. Physiology & Behavior 18, 283 (1977). 

[7] Y. Ashkenazy etal, Phys. Rev. Lett. 86, 1900 (2001). 

[8] P Ch. Ivanov etal. Nature 399,461 (1999). 

[9] M. Golubitsky etal. Nature 401, 693 (1999). 
[10] Y. Q. Chen etal, Phys. Rev. Lett. 79, 4501 (1997). 
[11] C. C. Chow etal. Hum. Movement Sci. 18, 725 (1999). 



